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Abstract 

In the framework of the recently proposed QCD based parton model for the 
cumulative phenomena in the interactions with nuclei the dependence of the cumu- 
lative pion production rates on the transverse momentum K± is studied. The mean 
value of K± is found to grow with x in the cumulative region. The obtained results 
are in agreement with the available experimental data. 
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1 Introduction 



A few years ago we have proposed a quark-parton model of cumulative phenomena in the 
interactions with nuclei |l|, based on perturbative QCD calculations of the corresponding 
quark diagrams near the thresholds, at which other quarks ("donors") in the nuclear 
flucton transfer all their longitudinal momenta to the distinguished active quark and 
become soft. 

Consider the scattering of a hadronic projectile off a nucleus with the momentum P 
in the cm. system. At high energies the momentum K of the produced pion belongs 
to the cumulative region if > Pz/A. As a reasonable first approximation, we treat 
the nucleus as a collection of = 3A valence quarks, which, on the average, carry each 
longitudinal momentum xqPz/A with xq = 1/3. In our approach the cumulative pion 
production proceeds in two steps. First a valence quark with a scaling variable x > 1 is 
created. Afterwards it decays into the observed hadron with its scaling variable x close 
to the initial cumulative quark's one. This second step is described by the well-known 
quark fragmentation functions and will not be discussed here. 

The produced cumulative ("active") quark acquires the momentum much greater than 
XqPz/A only if this quark has interacted by means of gluon exchanges with other p quarks 
of flucton ("donors") and has taken some of their longitudinal momenta (see Fig.l). If 
this active quark accumulates all longitudinal momentum of these p quarks then Kz = 
{p + 1)xqPz/A and the donors become soft. It is well-known that interactions which 
make the longitudinal momentum of one of the quark equal to zero may be treated by 
perturbation theory 0]. This allows to calculate the part of Fig. 1 responsible for the 
creation of a cumulative quark explicitly. This was done in Q , to which papers we refer 
the reader for all the details. As a result we were able to explain the exponential fall-off 
of the production rate in the cumulative region. 

Since with the rise of x the active quark has to interact with a greater number of 
donors, one expects that its average transverse momentum also grows with x. Roughly 
one expects that {K'j} is proportional to the number of interactions, that is, to x. In 
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H, H this point was not studied: we have hmited ourselves with the inclusive cross- 
section integrated over the transverse momenta, which lead to some simplifications. The 
aim of the present paper is to find the pion production rate dependence on the transverse 
momentum and the mean value of the latter as a function of x in the cumulative region. 
This dependence and also the magnitude of {K]_) have been studied experimentally. The 
comparison of our predictions with the data allows to obtain further support for our model 
and fix one of the two its parameters (the infrared cutoff). 



2 The dependence 

Repeating the calculations of the diagram in Fig.l described in |l|, |^ but not limiting 
ourselves with the inclusive cross-section integrated over the transverse momentum, we 
readily find that all dependence upon the transverse momentum K_i_ of the produced 
particle is concentrated in a factor: 

J{K^) = J pA{r^\r^^)G{ci, ...,Cp)f[X{c, - r)\{c, ^ (1) 
p+i p+i 

Here pA is the (translationally invariant) quark density matrix of the nucleus: 

- f _ ^ 

pA[ri\ri) = I ip^A{ri,r.ra)i^lAii^i,rm) H ^^^"^ (2) 

m=p+2 



where i{j±a is the transverse part of the nuclear quark wave function. The propgation of 
soft donor quarks is decribed by 

A(c) = !^ (3) 

ZTT 

where m is the constituent quark mass and Kq is the modified Bessel function (the Mac- 
Donald function). The interaction with the projectile contributes a factor 

p 

n 

i=l 

where (Jqq{c) is the quark-quark cross-section at a given value of impact parameter c and 



G{ci, ..■,Cp) = J Yl (^qq{ci - bi)r]Hibi, bp)(fbi (4) 



VH{bi,...,bp) = Y: / \^±H{h)\'5^'\^j:h) d\^,...d%L (5) 



is a multiparton distribution in the projectile, expressed via the transverse part of its 
partonic wave function ip^u . If we integrate J{Ki_) over we come back to our old 
result (Eq. (33) in 



P+i p+i 



i=l 



If one assumes factorization of the multiparton distribution rinibi, bp) then 
G{ci, ...,Cp) also factorizes: 

G{c,,...,Cp) = f[Go{c,) (6) 

i=l 

Following PI we use the quasi-eikonal approximation for rjH- 



i=l 



where ^ is the quasi-eikonal diffraction factor, u^j is the mean number of partons in the 
projectile hadron and the single parton distribution rj^ib) is normalized to unity. In a 
Gaussian approximation for a{c) and r]H{b) we find: 

where cjgg is the total quark-quark cross-section, r^^ = fQ+r'jj, Tq and rij are the widths 
of ct(c) and rjuib) respectively. 

With the factorised G(ci, ...,Cp) (|]) we have 

J{K^)= j pAiO, -,017 - r, ...,r - r)f{r,r)e'^^-''^'^^d\d'^r 
P+i p+i 

where 

j{r,r) = / d'^cGo{c)X{c — r)X{c — f) 



We also have used the translational invariance of the p-matrix. Note that near the real 
threshold we have no spectators and 

Pa(0^ I r-r,...,r-r) = p^(0^ | 0^) 

p+i p+i p+i p+i 
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In any case large i^^ corresponds to small r — r so we factor out of the integral in zero 
point. In the rest integral we pass to the variables 

B = , = r — r 

2 ' 

and shift the integration variable c, then 

JiK^) = Pa(0^ I 0^) I fiB, b)e''''^d%d^B (7) 
P+i p+i 

where 

jiB, b) = J GoiB + c)A(^ - c)A(^ + c)d'c (8) 

3 The calculation of 

Now we would like to find the width of the distribution on K± as a function of p or what 
is the same of the cumulative number x = {p+l)/3. From the mathematical point of view 
it is simpler to calculate the mean squared width of the distribution {K']). Unfortunately 
in our case this quantity is logarithmically divergent at large K^. This divergency results 
from the behavior of h) at small h. This behavior is determined by the behavior of 
the A(6) = -K'o("^|&|)/(2vr) (j^), which has a logarithmical singularity at \b\ = 0. Smooth 
Go{B + c) does not affect this behavior. 
For this reason we shall rather calculate 

= i-/ f{B,b)\K^\e^'^-dVB^ (9) 

where Jn is the same integral as in the numerator but without \K±\. Presenting \K±\ as 
K\/\K_[_\ and K\ as the Laplacian applied to the exponent we find 



(\K^\) = —r I f{B,b)A,e'''''-d%d'B- 



>N 

Twice integrating by parts and using the formula 
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l^±l \b\ 



we find 

Now we again integrate by parts once to find 

where = h/\h\. Tliis leads to our final formula 

(I^J) = / d'B'^^f-\BMn,V,)j{BM (10) 

where i{B, h) is given by (||), \{h) is given by (H) and 

7^ = 1 d^Bf{B,h = Q) 

4 Approximations 

To simplify numerical calculations we make some additional approximations, which are 
not very essential but are well supported by the comparison with exact calulations at a 
few sample points. 

As follows from the the asymptotics of Kq{z) at large z 




the width of A (6) (j^) is of the order m~^. The function Go is smooth in the vicinity of 
the origin and its width r^H = \/i^q + i^h is substancially larger than the width of A. For 
this reason we factor Gq{B + c) out of the integral (||) over c at the point B: 

j{B, b) = Go(5)A(&), A(6) = I A(c)A(c - h)£c = -^K^{m\h\) (11) 

J Arcm 

Then we find that the integrals over B and b decouple 

J{K±) = pa{0, I 0, 0) J GI{B)£B J k^{h)d^^^d% (12) 
p+i p+i 

In this approximation we find that (l-ft'xl) depends only on one parameter - the con- 
stituent quark mass m, which in our approach plays the role of an infrared cutoff: 

i"00 

{\K^\)=pm dzKo{z){zK,{z)r-^ (13) 
Jo 
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This allows to relate m directly to the experimental data on the transverse momentum 
dependence. 

5 Comparison with the data and discussion 

The integral in (|l^) can be easy calculated numerically. For values of p = 1, 12 it is 
very well approximated by a power dependence (see Fig.2), so that we obtain 

(|irx|)/m = 1.594 /-^^s (^4^ 

As we observe, the rise of (l-ft'xl) turns out to be even faster than expected on naive 
physical grounds mentioned in the Introduction (~ y/p)- The resulting plots for (IfTxl)^ 
as a function of the cumulative number a; = (p + l)/3 at different values of parameter m 
are shown in Fig. 3 together with avaiable experimental data from |^ on {Kj_) for pion 
production obtained in experiments with 10 GeV protons and p, ^ with 8.94 GeV 

protons. 

Note that earlier publications of the first group |^ reported a much stronger increase 
of {Kj_) with X, up to value 2 {GeV/cY at a; = 3 for pion production. In our approach such 
an increase would require the quark mass to be as high as m ~ 225MeV. In a more recent 
publication the rise of {KjJ is substancially weaker (it corresponds to m ~ 175MeV 
in our approach). The authors of ||^ explain this by new experimental data obtained and 
by a cutoff K±max introduced in calculations of {Kj_) in ||^. The introduction of this 
cutoff considerably (approximatly two times) decreases the experimental value of {K'j_) at 
a; = 3. In our opinion this is a confirmation that the cumulative pion production rate only 
weakly decreases with K± in the cumulative region so that the the integral over Kj_ which 
enters the definition of {KjJ is weakly convergent or even divergent, as in our approach. 
Undoubtedly presentation of the experimental data in terms of the mean value (IK^I)"^, 
rather than {Kj) should reduce the dependence on the cutoff K±rnax and make the results 
more informative. 

One of the ideas behind the investigations of the cumulative phenomena is that they 
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may be a manifestation of a cold quark-gluon plasma formed when several nucleons overlap 
in the nuclear matter. In [Q we pointed out that our model does not correspond to this 
picture. It implies coherent interactions of the active quark with donors and, as a result, 
strong correlations between the longitudinal and transverse motion. Predictions for the 
dependence of (l-ft'xl) on x are also different. From the cold quark-gluon plasma model 
one expects (l-ft'xl) to behave as x^^^, since the Fermi momentum of the quarks inside the 
overlap volume is proportional to the cubic root of the quark density. Our model predicts 
a much faster increase, with a power twice larger. The experimental data seem to support 
our predictions. 
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Figure captions 

Fig. 1 The diagram for the production of a cumulative quark with the momentum K in 
the scattering of a projectile hadron with the momentum if off a nucleus A with 
the momentum P. Dashed and chain lines show gluon and pomeron exchanges, 
respectively. 

Fig. 2 The {\K^\)/m as a function of p. The points are the results of calculations on 
(p!3|). The line is the best power fit. 

Fig. 3 The {\K^\Y as a function of the cumulative number x = (p + l)/3. The lines are 
the results of calculations on ([T3|) at different values of parameter m. The points (•) 
are the experimental data from ^ on {K"]) for pion production with a cutoff (see 
the text) obtained in experiments ^-[^ on a bombardment of nuclei by 10 GeV 
and 9 GeV protons. The points (o) are the data from the earlier publications of the 
group [D, without a cutoff. 
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